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Abstract 

The observed results for the azimuthal single spin asymmetries (SSAs) of the proton, measured in 
the semi-inclusive deep inelastic scattering (SIDIS), can be explained by the final-state interaction 
(FSI) from the gluon exchange between the outgoing quark and the target spectator system. SSAs 
require a phase difference between two amplitudes coupling the target with opposite spins to the 
same final state. We have used the model of light front wave functions (LFWFs) consisting of 
a spin-^ system as a composite of a spin-^ fermion and a spin-1 vector boson to estimate the 
SSAs. The implications of such a model have been investigated in detail by considering different 
coupling constants. The FSIs also produce a complex phase which can be included in the LFWFs 
to calculate the Sivers and Boer-Mulders distribution functions of the nucleon. 
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I. INTRODUCTION 


In recent years, there has been lot of interest to investigate the single spin asymmetries 
(SSAs) PiE] which are considered the most interesting phenomenon to understand Quantum 
Chromodynamics (QCD) from the basic principles. When the target nucleon is transversely 
polarized to an incoming beam, many particles (hadrons) are produced in collision which 
tend to show right-left asymmetry in their distribution relative to nucleon spin direction. 
One of the participating particles in scattering process carries a polarization and if the 
scattering cross section depends on the direction of this polarization, a SSA is obtained. 
One of the remarkable aspect about the SSAs is that they persist at high energy and show a 
very stable pattern. A strong correlation between the target proton spin Sp and the plane of 
the produced pion and virtual photon has been observed in the semi-inclusive deep inelastic 
scattering (SIDIS) experiments by the HERMES and SMC collaborations |3H3]. COMPASS, 
a hxed target experiment at CERN |6Hin!, presents the result for the SSA of charged hadrons 
produced in deep inelastic scattering of muons on a transversely polarized target. STAR 
collaboration [nHl3] shows that the SSA for pp —>■ tiX process at center of mass energies 
of 20 and 200 GeV may arise from the Sivers effect [1], Collins effect [H] or a combination 
of both. A large azimuthal SSA was observed in the hadronic reactions like pp^ —)■ tiX 
[15] and in pp —)■ A^X pTO] where the antiproton and hyperon respectively are polarized 
normal to the production plane. The issue regarding SSAs becomes all the more difficult to 
understand when it is realized that the sign of the SIDIS experimental and Drell-Yan (DY) 
experimental im measurements are different. 

The fact that the understanding of SSAs theoretically depends simultaneously on several 
aspects of hadron structure, for example, interference between complex phases, orbital an¬ 
gular momentum and hnal state interactions IFSIH [171120] . makes the study of SSAs even 
more interesting. It was shown in Ref. [21] that a relation exists between the generalized 
parton distributions (GPDs) and SSAs (SSAs= ESI X GPDs). In particular, the transverse 
distortion of the distribution of quarks in a transversely polarized target nucleon is related 
to the spin flip GPD E{x,0,t). It was shown that FSIs from gluon exchange between the 
outgoing quark and the target lead to SSAs [IH| which essentially require a complex phase 
difference between the two amplitudes which couple the proton target with = ±| to the 
hnal state. The SSAs have been studied for both the SIDIS and DY process by using a QGD 
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motivated quark-diquark model of a spin-| proton of mass M with spin-^ and spin-0 con¬ 
stituents of mass m and A respectively |22]- SSA for DY process was explained with the help 
of initial state interactions (ISI) [171 [23] • Collins [Tl| pointed out important consequences of 
SSAs in deep inelastic scattering (DIS) by incorporating the struck quark into Wilson line 
path order exponents. The theoretical description of the single spin phenomena in QCD is 
still a big challenge. 

A comprehensive picture of the nucleon can been obtained by considering the transverse 
momentum dependent parton distributions (TMDs). TMDs are the T-odd parton distribu¬ 
tions which give rise to SSAs and describe various hard exclusive reactions of the nucleon 
[211125]. The FSIs can also act as a source of the Sivers and Boer-Mulders distribution func¬ 
tions [261 - 128] . Sivers distribution function proposed by Sivers in 1990 [1], describes the 
distribution of unpolarized quarks inside the nucleons which are transversely polarized in the 
opposite direction. Due to PT invariance Sivers asymmetry was supposed to be prohibited 
/it = -/it = 0 [I3|. However, Brodsky, Hwang and Schmidt proved that /ir ^ 0 [TEj. On 
the other hand, Boer-Mulders function (h/) describes the distribution of polarized quarks 
inside the unpolarized nucleons. Extensive amount of work has been done in the light cone 
quark models [29113T] as well as diquark models which include the inclusive pion and kaon 
production in DIS [12] . Sivers distribution function has also been discussed in the MIT bag 
model with the help of interference of S and P wave components in the presence of FSIs 
[32]. Further, model independent and dependent relations between the GPDs and TMDs 
have also been developed [33] . The derivatives of chiral odd GPDs in the impact parameter 
space can be related with and h/ as well as with the spin densities of the nucleon [33]. 
The SSA and the Sivers distribution function can be related on the basis of the interference 
terms coming from the FSI. Recently, Sivers and Boer-Mulders distribution function have 
been studied on the lattice to non-local operators [35]. We would like to emphasize that 
the signs of and h/ are of great interest as different models predict different signs for 
them. Therefore, it becomes interesting to study the Sivers and Boer-Mulders distribution 
function in the context of SSAs and FSIs. 

The set of light front wave functions (LFWFs) provide a frame-independent, quantum- 
mechanical description of hadrons at the amplitude level which are capable of encoding 
multi-quark and gluon momentum, helicity and flavor correlations in the form of universal 
process independent hadron wave functions. One can also construct the invariant mass 


3 


operator Hlc = P^P~ — P± and light cone time operator P~ = in the light cone 

gange from the QCD Lagrangian [SBtiHS] . Since the essential physics that allows the SSA 
is the existence of two different light front angular momentum components in the proton 
Lz = 0, Lz = PI with different FSIs, therefore it becomes interesting to calculate the SSA 
in the SIDIS process 7 *p —)■ q{qq)i where a virtual photon is scattered on a transversely 
polarized proton having transverse spin and is induced by the FSIs. In this model, the spin- 
I proton of mass M is a composite of spin-| fermion of mass m and spin-1 vector boson of 
mass A and it lies in the framework of QED (SHHll] • This model has been successfully used 
as a template for calculating the Schwinger anomalous magnetic moment, understanding the 
structure of relativistic composite systems and their matrix elements in hadronic physics, 
explaining the hadronic structure and providing a good representation of 3-D structure 
of the nucleon gni SB m]. Further, it also gives the general proof for the vanishing of 
the anomalous gravitomagnetic moment 5(0) [3B SB S3]- This model has been used to 
calculate the spin and orbital angular momentum of a composite relativistic system as well 
as the GPDs in impact parameter space [151 SB] . 

In the present work, we have calculated the single spin asymmetry (SSA) in the semi- 
inclusive deep inelastic scattering (SIDIS) process 7 *p —)■ q{qq)i. The amplitude for this 
process is computed at both the tree and one-loop level and the SSA has been calculated by 
taking interference between these amplitudes. Final-state interaction (FSI) from the gluon 
exchange between the outgoing quark and the target spectator system has been used to 
explain the observed results for the azimuthal SSAs of the proton. SSAs require a phase 
difference between two amplitudes coupling the target with opposite spins to the same hnal 
state. To estimate the SSAs, we have used the model of light front wave functions (LFWFs) 
consisting of a spin-| proton of mass M as a composite of a spin-| fermion of mass m and a 
spin-1 vector boson of mass A. The FSIs also produce a complex phase which can be included 
in the LFWFs to calculate the Sivers and Boer-Mulders distribution functions of the nucleon. 
Sivers function is given by the overlap of the wave functions having opposite proton spin 
states but the same quark spin state whereas the Boer-Mulders function has the same proton 
spin state but opposite quark spin state. The Sivers and Boer-Mulders distribution functions 
have been computed by inducing the spin-dependent complex phases to the LFWFs. The 
implications of such a model have been investigated in detail by considering hxed and running 
coupling constants. 
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FIG. 1. The (a) tree diagram and (b) one-loop for 7 *p —>■ q{qq)i. The interference amplitudes 
obtained from these two diagrams provides the SSA. 


II. LIGHT FRONT WAVE FUNCTIONS (LFWFS) 


The single spin asymmetry is calculated in the SIDIS process 7 *p —)■ q{qq)i where virtual 
photon is scattered with virtuality = —Q^ on a transversely polarized proton having 
transverse spin where a spin-| system acts as a composite of a spin-| fermion and a spin- 
1 vector boson. During the process, the virtual photon is absorbed by an active spin-^ 
quark as shown in Fig. [T] and the diquark is used to describe the spectator system. The 
tree level diagram (Fig. [^(a)) alone will not describe the observed asymmetry because the 
contribution is real towards the scattering amplitude. However, the interference between 
the tree level and the one-loop (Fig. (b)) amplitudes together will generate the desired 
SSA. Kinematics have been developed in the Drell-Yan-West frame which is boosted so that 
g"*" = 0 and is collinear to proton [TH]. The longitudinal momentum fraction exchanged is 
dehned in terms of the momentum carried by the outgoing quark r and the momentum 
carried by the photon g as A = ^. We can fix r~ and q~ from the momentum conservation 
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and the on-shell conditions for the nncleon, qnark and diqnark and we have 

M2 + A2 


r = 


p+ (1 — A)p+ ’ 

+ 2g_L ■ f± 

a ~ -. 

Ap+ 

The other kinematical details can be snmmarized as follows 

^ =(^ 0^), 

((f_L -I- rj _)2 m 2 - 1 - A 2 


( 1 ) 


= (O, 


, m2 
^ 


Ap+ p+ (1 — A)p+ 
r-^ + A2 _ 


, , 


p+ (1-A)p+’ 

We consider here a spin-^ system which acts a composite of spin-| fermion and spin-1 vector 
boson [ElEHlIlSlSe]. If we consider x as the fraction of momentnm transferred and fc_L as 
the transverse component of the momentnm, the two-particle Fock state for a fermion with 
= -t-|, having fonr possible spin combinations, can be expressed as 

d^k I da; 


two particle 




\Jx{l — a;)167r3 


(3) 


X 




+ - + 1; a:P+ , _^(a;, A;±) 


+ - - 1; xP+ , 


+'<p\^^{x,k^) 


where 


^ -hi; xP+, _^(a;, fcj_) 


V'ti , fcx) = -\f2 (m - —] (p, 

2 \ X y 

i^^_i_Ax,k±) = 0. 

2 

Similarly, for ^ we have fonr possible combinations 


- — 1; xP^ , k± 
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^Loparticle(^^,^± = 0j 


d2fc I da; 


A/a;(l — a;)167r3 


(4) 


( 5 ) 


X 




+ - + 1; xP+ , A;± ) + _^(a;, k±) 


+ - - 1; xP+ , fcx 


+^pt^^^{x,k^) 


-- -hi; xP+, A;± )+^/;_i _^(x, fcx) 


■ - — 1; xP~^ , fcj_ 
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where 


= 0, 


ki_) = -\f2 (m - —] if, 

2 \ X / 


fc±) = -V2- 


-k^ + ik"^ 


(i-i) *’• 

i’*'t (i.fci) ^ -V2 * ^ 'f 

~2-^ x[l — x) 


and 


(^(x, fcj_) = 


\/r 


X 


M2- 


/c^+m^ fc^+A^ 

X 1—x 


( 6 ) 


(7) 


Here the formalism has been generalized by assigning a mass M to external fermions in 
the scattering process bnt a different mass m to the internal fermion line and a mass A 
to the internal boson line. The charge of the fermion is taken as e. It may be important 
to mention here that the nnmerators in the wave fnnctions are characteristic of the orbital 
angnlar momentnm and hold for both pertnrbative and non-pertnrbative conplings. 

III. TREE AND ONE LOOP AMPLITUDES 


The SSA can be calcnlated from the interference of the amplitndes at tree and one-loop 
levels as shown in Fig. The tree level amplitnde, calcnlated from Fig. (a), receives 
contribntion from the following amplitndes with an 0{as) FSI 

= = (-r‘ + *r^)C 

St = -l) = -V 2 


A^^^%i^^l,st = +l) = -V2 

= -1) = -V2 


A(1 

-A)(ri 

-M2A(1-A)+m2(l- 

- A) + A2A)’ 



{A -\- iA)C 



(1- 

A)(ri- 

M2A(1-A)+m2 

’(1- 

A) + A2A)’ 

:(^M 

m\ 

U 



a) {r‘j_ 

-M2A(1-A) + 

m2(l 

- A) + A2A) 

:(^M 

m\ 

C 



a) {r‘j_ 

-M2A(1-A) + 

m2(l 

- A) + A2A) 



{—A + iA)C 



A(1 

-A)(ri 

— M2A(1 — A) -|- m2(l - 

- A) + A2A)’ 



{A + iA)C 



(1- 

A)(ri- 

M2A(1-A)+m2 

’(1- 

A) + A2A)’ 


where 


C = -e eiP+yA2 A (1 - A). 


( 8 ) 


(9) 
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We have taken here the electric charges of q and {qq)i to be ei and 62 - The state ft (If) 
denote the two spin states of the proton whereas for the spin projection of the 

spin-^ constituent = ±| we label t and j,. 

The one-loop amplitudes can be expressed as 


A°n^-i°°P{I) = iee\e2 


d^k M{I) 


— m? + ie){{k + — m? -|- ie) 

1 

{{k — rY — -I- ie){{k — pY — -|- ie) ’ 


X 


( 10 ) 


where Xg is the mass of the gauge boson. In the standard leading order perturbative QCD 
(pQCD) calculations |17], hnal state multiple interactions are considered whereas in the 
present calculation the spin asymmetry is obtained by considering a single hnal state inter¬ 
action. The spin asymmetry in pQCD is obtained from the non-vanishing transverse momen¬ 
tum dependent cross section and in our calculations we consider the process 7 *p —)■ q{qq)i 
where the detected particle is identical to the quark and the asymmetry for this detected 
hadron can be obtained by convoluting the jet asymmetry result with the realistic fragmen¬ 
tation function. A similarity in both the approaches is that a required phase is essential to 
generate the asymmetry and this phase comes from a pole in the propagator. 


The numerators A/'(/) in Eq. (10) can be obtained from the Feynman diagram (Fig. 
(b)) as follows 


(si = +1)) = 2P^VA x( - 

(Sfo = -1)) = 2P+yA x(^-V2 

{si = +1)) = 2P+^/A (m - 

Af(^^t, {si = -1)) = 2P+^/A (m - 

Af(^^t, {si = +1)) = 2P+VA x( - V2 

{si = -1)) = 2P+yA x( - V2 (11) 


where q = 


_ 2 !_ 

AP+ 


with A = ^ being the quark light cone fraction. The denominators in 


p+ 










Eq. (10) can be rewritten as 


— ie 


— m? + ie = k'^ (k - 

V fc+ 

(k + qf - + « = (k* + ,+)((i- + ,-) - - “ 


(fc — r)^ — + ie = (fc‘ 
(A; — p)^ — + ie = {k~ 


)((fc -r ) 

(A;--P-) + 


xP+ 

Ag + (A;± - r±f - ie 
(a; - A)P+ 

A2 + A;2 - ie 


( 12 ) 


(1 — a;)P+ / ’ 

whereas the integral over the momentum space in light front held theory can be expressed 
as 

f d'^k = 


d‘^k I dk^dk 


The one-loop amplitude now becomes 


jone-loop^j-^ = -iee^es 


J 2(27r)4 

d^k^ r P+Af{I)dx 

2{27iy J {P~^)‘^x‘^{x — A){1 — x) 

dk~ 


(13) 


X 


ie 

xP+ 


)((*:- +r) - 


X 


1 


{k~ — r~) 




(^- - P-) + "PAP 


(14) 


(a:-A)P+ 

We would like to mention here that there is a restriction on the combination of propagators 
that can go on-shell simultaneously and the cuts have to be imposed by hand. These 
kinematical constraints have been discussed in detail in Ref. [22]. On performing the 
integration over k~, we get 


^one-ioop^j^ = -iee?e 2 ( 27 ri) 


d^kx 

2(27r) 


P^ 


U{I)d. 


X 


{P+Yx‘^{x — A)(l — x) 
1 


X 


p- - 


—ie —ie 


(1—x)P+ 


a:P+ 


_ A^+fc^—ie m2+(fcx+qx)^— 


- 7T3 


(1—a;)P+ 


xP+ 


X 


P- 


—ie A2 + (Aix— rj^)2—i 


(16) 


(1-3;)P+ (3;-A)P+ 

The imaginary part of the above equation gives rise to the imaginary phase which is essential 
for the SSAs. This imaginary part evolves from the real intermediate propagator state 
allowed just before re-scattering and can be dehned in terms of the initial state (PPit) and 
intermediate state {PAterm) energy denominators as 

1 


Im 


P — P 7 f 

init interm ' 


=m 


init 


interm) 


(16) 
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The imaginary part from the propagator term from Eq. (15) gives 


X^+k^—ie rn'^+{k±-\-q±)'^—ie 
xP+ 


P~ + q~ 


{l—x)P+ 


= —iird yP + q 

^2p+ 


+ k‘]_ rrP + {k± + qx)" 


{l-x)P+ 


xP+ 


= —ITT- 




6{x-A-S), 


(17) 


where 


5 = 2A 


q± ■ {k± - fj_) 




(18) 


In terms of <5 the exchanged momentum 6 P~^ is very small. This leads to the light cone 
energy denominator being dominated by term. 

Using Eqs. (0, 0 ,and ([TT]), one-loop amplitudes obtained can be written as 


j^one-loop^^^^^ (Sg = +1)) 

j^one-loop^^^^^ (S^ = -1)) 

A°ne-l°°P(^^l^ (sl = +1)) 

A°n^-l°°P{i^^t, {4 = -1)) 

{sl = -1)) 

^one-loop (Sg = +1)) 


-iee^e2 


iee^e2 


iee^e2 


iee^e2 


—166162 


iee^e2 


r d^kx2P+x/A A (1 - A)V2{k^ - ik^ 

I 167r2A(l - A)LiL2 ^ 

d^kx2P+x/A A (1 - A)V2{k^ + iP) 

167r2(l - A)LiL2 ’ 

d^k±2P+x/A A (1 - A)V2{M - f) 
L1L2 

d^kx2P+^ A (1 - A)V2{M - f) 
r d‘^kx2P+y/A A (1 - A)V2{k^ - ik^ 

I 167r2(l - A)LiL2 ' 

d^kx2P+VA A (1 - A)V2{k^ + ik4 
167r2A(l - A)LiL2 ’ 


(19) 


where 


Li = kl- M^A{1 - A) + m^{l - A) + A^A, 

L2 = \l + {kx-Txf. (20) 


These can now be integrated over the transverse momentum kx using Feynman parametriza- 
tion. Using Eqs. (Isj) and (19), the interference between the tree level and one-loop ampli- 
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tudes can be calculated and we get 

(-Sfe = + 1 )) = 


-\/ 2 (—+ ir'^] 

A(l-A) 




= -!)) = 


1 - A 


Ibvr^ 

. 6162 

IGtt^ 
. 6162 


(-Sfe — + 1 )) — —+ 

(4 = -!)) = 5)c('‘ + -16^2 

(,. = + 1 ), = zAtA^cffc + 


1 - A 


IGtt^ 


(si = + 1 )) = 


IGvr^ 


( 21 ) 


A(l-A) 

where C has already been dehned in Eq. 0 and 
/■' 1 

=i 

^ ^ - M2A(1 - A) + m2(l - A) + A2A ’ 

D = a{l-a) + a \l- M2(1 - a)A(l - A) + 771^(1 - a)(l - A) + A2(1 - a)A.(22) 


IV. SINGLE SPIN ASYMMETRY (SSA) 

We now present the results for the SSA. Following Refs. [EiEg we have used the gauge 
particle as an infrared regulator and have taken the mass of gauge boson as \g = 0. The 
calculations have been performed in the region A < a; < 1. The production plane has been 
taken as the z — x plane which can be dehned from the virtual photon and the produced 
hadron. The asymmetry is produced in the y direction and can be expressed as 

'^y ^ Q )-t); (sfe = +1) — {si = +1)) + 

7(A('fl —{si = +1) — A(ft'—>- 4 ,)*; {si = +1)) + 

{si = - 1 ) - {si = - 1 )) + 

7(A('fl—)-4,)*A(Jj.—)-4,); {si = —1) — A('fl—>'4,)A(1].—>- 4 ,)*; (s^ = —1)) j, (23) 

where 

C {si = +1))P + |A(ff'—)-t, {si = —1))P + )-4-, {si = +1))P + 

)-t, {si = —1))|^ + )■!, {si = +1))|^ + {si = —1))|^. (24) 
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Substituting the amplitudes from Eq. (21), we get azimuthal spin asymmetry transverse to 
production plane as follows 


eies r\MA - m)A(l - A)(r^ - - A) + rn^^il - A) + A^A) 

~ (r^(l + A2) + (MA - m)2(l - A)2) ^ 

1 - M2A(1 - A) + m\l - A) + A^A) 

"" -M2A(1 - A)+m2(l - A) + A2A ‘ ^ ’ 

Here the factor reflects that single spin asymmetry is proportional to the term 


Sp ■ q X f. 

The ESI from gluon exchange in the present case has the strength —)■ where 

Cp is taken to be | so that the results match with QCD. Momentum transfer carried by the 
gluon Exes the scale for the strong coupling constant ai which has a hxed value of cti = 0.3 as 
obtained in MS scheme [IH]. In the present work the masses of the nucleon, quark, diquark 
and gluon have respectively be taken as M = 0.94 GeV, m = 0.3 GeV, A = 0.8 GeV and 
Ag = 0 GeV [iH]. For the sake of comparison as well as a deeper understanding of the role 
of coupling constant we have also calculated the asymmetry by replacing the hxed coupling 
constant (ai) by running coupling constants where higher order contributions, particularly 
the closed quark loops, can be taken into account. We take an analytic one-loop running 
coupling constant as proposed by the Shirkov and Sovlovstov IIHHSI] 


«2(Q') 


dvr 1 

A HQVI^IcdY 


( 26 ) 


with = 11 — “^Nf {Nf being the number of active quarks). The use of such running 
coupling constant was found to be essential for explaining the reproduction of the light 
pseudoscalar meson spectrum. However, for Aqqd = an nonphysical singularity existed 
which contradicted some analytical properties and had to be modihed in the infrared region. 
Therefore, the above said coupling constant was modihed [HHST] such that it remains finite 
at Aqijd = and is given by 


A (ta(QVA|oD) ^ 

In order to check the behaviour and dependence of asymmetries on light cone momentum 
fraction A and momentum carried by the outgoing quark r±, in Fig. [^we have presented 
the model prediction for the transverse azimuthal spin asymmetry Vy as a function A and 
r± for the hxed coupling constant ai = 0.3 and the running coupling constants, 0^2 and 0 ^ 3 . 
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The longitudinal azimuthal spin asymmetry can be determined from the transverse az¬ 
imuthal spin asymmetry Vy using the relation 


A sin 0 
^UL 


KA 


sin 0 
UT 5 


( 28 ) 


where A^^ is HERMES transverse asymmetry where the target polarization is transverse 
to the incident lepton direction, is the HERMES longitudinal asymmetry where the 

target nucleon is polarized along the incident lepton direction |3]. The kinematic factor K 
is defined in terms of the virtual photon fractional energy y as 


K 


2Mx ll — y 

Eiab V y 


(29) 


We have taken y = 0.5 and Eiab = 27.6 GeV [52]. The longitudinal azimuthal spin asym¬ 
metry KPy is plotted in Fig. [^as a function A and r± for all the three coupling constants 
mentioned above. 


In Figs. (a) and (a) we have shown the variation of transverse and longitudinal 
asymmetries (Py and KVy) with A at a hxed value of momentum carried by the outgoing 
quark r_i_ = 0.5 GeV. Interestingly, the transverse as well as the longitudinal asymmetries 
change sign at A ~ 0.3 as the value of the A increases. This may be possibly due to the 
qq pairs being readily produced at this momentum fraction which give rise to sea quarks 
affecting the sign of asymmetry at this point. The sea quarks seem to play an important 
role in carrying the momentum fraction. Future experiments like E906/SeaQuest Drell-Yan 
experiment [53l |5l| are ready to study the sign change of asymmetries at the momentum 
fraction A ~ 0.3. Further, in Figs. (b) and (b) we have hxed the momentum fraction 
A = 0.15 GeV and plotted Vy and KVy as a function of quark transverse momentum r^. 
Since r±_ is the transverse momentum of the outgoing quark relative to photon direction, the 


asymmetries do not change sign in these cases. It is also clear from the Eq. (25) that Vy 
decreases as cxi A 0-/n r^. The values of and in the above calculations are taken to 
be 0.25 GeV and 0.10 GeV respectively. 
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FIG. 2. Transverse asymmetry Vy as a function of A and r_L of the proton obtained by gluon 
exchange in FSI. 




FIG. 3. Longitudinal asymmetry KVy as a function of A and r±_ of the proton obtained by gluon 
exchange in FSI. 

V. SIVERS AND BOER-MULDERS DISTRIBUTION FUNCTIONS 

Quark transverse momentum dependent parton distributions are defined through the 
correlation function as follows [55H57] 

t’nU.L.s) = |L|^e“'.(F,s|vyo)rwv-(!/)|P.s>|,*.„. (30) 


Here P is the nucleon momentum, x and are the fractional longitudinal and the transverse 
momentum of the quark, S is the covariant spin vector having components , 

S~ = and The gauge link (Wilson line) W assures color gauge invariance of 

the correlators and F are the usual gamma matrices. By projecting the correlators onto the 
full basis of 7 matrices one can obtain the corresponding distribution functions. At leading 
twist the proton spin structure can be described by 8 TMDs. Using the gamma matrices 
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7’*', 7'''75 and i(T-^’'' 75 , we obtain 




^ frOp ci 

5'Lfi'iL(a;, fcj_) + ^a), 

^ \U ^ , 

+ -^h^L^x.k^) + S^^\hiT{x,k^) + 

2 k{kx ■ s^ - 

iM 2 



+ 


(31) 


The LFWF representation of the Sivers {f^{x, k±)) and Boer-Mnlders {hj;{x, k±)) distribn- 
tion fnnctions can be dednced from the proton eigensolntion |pp) on the eigenstates {|n)} 
of the free Hamiltonian giving the light cone Fock expansion as 


wj’*T±))=En 

n,Aj 2=1 


dxi (^k±i 
Ibvr^ 


Ibvr^^ 



^( 2 ) 


X'lpniXi, fc_Li, \i)\n-,XiP^,XiP±_ + /c_Li, A*). 



(32) 


The relative momentnm coordinates of the LFWFs are the momentnm fractions x* = ^ and 
the transverse momenta k±_i of partons. The light cone spin projections along the direction 
of qnantization are expressed as Aj whereas the physical transverse momenta of partons are 
Pxi = XiP± + k±i [58] . 

The distribntion fnnctions fi{x,k±) and f^{x,k±) can be dehned throngh the matrix 
elements of the bilinear vector cnrrent as follows 


dy d?y 


IGtt^ 


J- ^ixP+y -ifcx-P 




2F+ L 


y+=0 

k\ 


h(x,hi) U{P,Si) 7+ U(P,Si) + f^-AxXi) g U(P,Si) o'* U(P,Si) 


(33) 


where 

^f/(F, 5^)a'+f/(P, S^) = with 6^2 ^ _pi ^ ^ _ ^ 34 ) 

For the present work, we consider the transverse spin S± = (Pi,*?!) = (0,1). The proton 
state can be dehned as a combination of two states as ^ind we have 

r dy d y± (-P; t I ~ i{P,l I p/'g'i l~^’ 

J 167r3 ^2 ^ 72 P=o 

= fi{x, ki_) - ^ fiT^x, fc±). (35) 
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We can now write 


h{^,ki)= 2 ^ 


(^’,tl^+(9)l^’,t) + {^’4l^+(9)l^’4) 


/lifeL) = ^ [{p.t 4+(!/)|/’4) - (^4 4+(!/)l-P.t>' 


y+=0 


y+=0 

(36) 


with J~^{y) = 'ijj{0)'y~^'ip{y). 

The Boer-Mnlders distribntion fnnction can be dehned throngh the matrix elements of 
the bilinear tensor cnrrent as follows 


dy d'^y 


Ibvr^ 


± ^ixP+y -ikj_-yj_ 


{P,S^\m a^^^l^{y)\P,S^) 


y+=0 


2P+ 


hi(x,k^)f^U(P,S^)7* U(P,S^_) 


(37) 


leading to 


k. 

M 




y+=0 

(38) 


with = '^(0)(T*+'^(|/). 

The light cone representation of nnpolarized qnark distribntion and Sivers fnnction can 


be written from Eq. (36) and is expressed as 


fi{x,k^) = C' - 


1 r 


k±i, Xi) k±i, Xi) + k±i, Xi) V’(n)(a^i, k±i, Xi) , (39) 


,t *1 


k] 


M 

where 


^ k±) = ^ ^In) k±i, Xi) k±i, Xi) - (xi, k±i, Xi) ijJ^){xi, k±i, Xi) 


i 


(40) 


^ =11 jti ( 1 - ] 5 ( 2 ) j ^ j _ 

n,Ai i=l \ j=l J \j=l J 

(41) 

On the other hand, the light cone representation of the Boer-Mnlders fnnction can be written 


from Eq. (38) and is expressed as 
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k] 


C 


* 


(^7 9 ^ \ "^1 4^; -^ZT^l) V^(n)(^^’ ^J-^? "^1 T: "^ 27 ^ 1 ) 


'^(n) ^-L 25 -^1 T? ^i^l) ^J-i; -^1 "^27^1) 


+ 


^( n ) ^ J - i : "^1 ' I '? -^ ZT ^ l ) ^ J - i ? "^1 T ? ^i^l) 


^(n) -^1 T? "^27^1) ^J-27 -^1 J ^7 -^27^1) 


( 42 ) 


It is clear from Eq. (40) that the Sivers function is the product of LFWFs having opposite 
proton spin states but the same quark spin state whereas for the Boer-Mulders function 


from Eq. (42) we note that LFWFs have same proton spin state and opposite quark spin 
state. 

For the calculations of Sivers and Boer-Mulders distribution functions in the model of 
spin-1 diquark model, we adopt the same treatment for the LFWFs as in the earlier sec¬ 
tion. We have induced spin-dependent complex phases to the wave functions which can be 
expressed as 


/t ( t \ -V2{-k^ + ik‘^) 


x(l — x) 

-V2{k^+ik^ 


1 -I- i-^ik‘]_ -I- B)l4 ) (p, 


1 -I- i-^{k‘j_ + B)!^ ) (p, 


fcx) = -\/2{m - ^) (l + + ^)^3) 

= 0 ) 

2 

= -V2(^M - ^) (l + + ^)^3) 

y^i-k^ + ik"^) 


(43) 


2 

/,i 


1 -I- i-^{kj_ -I- B)l4 ) (p, 


(1 -x) 

/r / ; 4 -V 2 {k'^+ik‘^) ( .£162,^2 , 7^4 r \ 

^_4_,(x, (1 + ^^{kx + B)h) 

= 0 , 


where 


(f = —e- 


X 


vT 


X 


P^+B 


(44) 

(45) 


with B = — M^x(l — x) -|- m^(l — x) -|- A^x and 


1.^1 ^da, 


(46) 
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FIG. 4. Results for Sivers and Boer-Mulders function for spin-1 diquark model. 


with 


D' = a{l — a) k‘]_ + a — M^(l — a)x{l — x) + m^(l — a)(l — x) + A^(l — a)x. (47) 


Using Eqs. (40), (42), (43) and (|44|), we can write 


fi{x,k±) = ( 




x‘^{i — xy 

—* / 777 / \ 

/i^(x, fcu) = 2M(m -j 

\ oc y 




(d+ 

x(l - x) Stt ^ ^ B 


m\ 1 + x 6162 


hf(x, ky = M(M - - ] ' \ (ki + BW^ln 

^ ^ ^ V x/x(l -x) Stt ^ ^ 


B 

k\ + B 
B ' 


(48) 


The azimuthal spin asymmetry Py can be expressed as a factorized form of the Sivers 
function as follows 


Py = 




IT 


M /. 


(49) 


In the present work, the azimuthal spin asymmetry Py and the Sivers function f^plx, k±) 


calculated respectively in Eqs. (25) and Eq. (48) are found to satisfy the above factoriza¬ 


tion. In the standard pQCD approach 03. the azimuthal spin asymmetry factorizes as the 
product of Sivers function, quark distribution function, hadronic cross-section and quark 
fragmentation function. The cross-section in this case factorizes as the product of quark 
distribution function and quark fragmentation function for SIDIS at large Q^. 

In Fig. I^we have presented the results for Sivers and Boer-Mulders distribution functions 
as a function of x for the fixed and running coupling constants ai, 02 and 03 . The value 
of k± is taken to be 0.1. We hud that both Sivers and Boer-Mulders distribution functions 
are overall negative for x ^ 0.3 but change signs for x ^ 0.3. Similar behaviour has 


already observed in the relation for asymmetry in Eq. (25). An overall negative sign for 
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the Sivers and Boer-Mulders distribution functions has already been obtained [HI |23l |58] • 
This occurrence of node in the present model is because of the term which clearly 

changes sign around x ~ 0.3 at the given value of masses in Section IV. The possibility of a 
node in the a:-dependence of the Sivers and Qiu-Sterman (QS) function T(a:, S'_l), related to 
each other by direct proportion = ^T{x, S±) |59l - l6T] . has also been discussed if 

the light cone integral is considered over the gluonic held strength [62j. Some other model 
calculations |63l |6l] have also shown the nodes for the up and down quark Sivers function. 

Calculations have also been done in the pQCD [65l |66] where Sivers asymmetry is shown 
to be related to the QS function by considering a point like and a dipolar form factor for the 
interaction between the nucleon, the quark and the spectator scalar and axial-vector diquark 
respectively. In Ref. iZlEH], the diquark model has been considered by taking the diquark 
as an axial-vector rather than a vector as considered in the present work. They have also 
considered the longitudinal polarization vector in addition to the transverse polarization 
vectors. The result obtained for Sivers and Boer-Mulders distribution functions have a sign 
difference between them but in our calculations no sign difference has been predicted. In 
addition to this, the overlap representation of Sivers and Boer-Mulders in Ref. [671168] is 
an assumption as it is not known a priori if the FSI operator can be isolated and the FSI 
operator is same for all functions and all types of diquarks. However, in our calculations 
the Sivers and Boer-Mulders distribution functions have been explicitly dehned from the 
overlap representation of wave functions. This can be substantiated by the results of several 
future experiments that are contemplating the possibility of performing the high precision 
measurements over a wide x region. 


VI. CONCLUSIONS 

We have calculated the single spin asymmetry (SSA) in the semi-inclusive deep inelastic 
scattering (SIDIS) process YP q{qq)i using the model of light front wave functions 
(LFWFs) consisting of a spin-| proton of mass M as a composite of a spin-| fermion of mass 
m and a spin-1 vector boson of mass A which lies in the framework of QED. The amplitudes 
for the SIDIS process are computed at both the tree and one-loop level. SSAs require a phase 
difference between two amplitudes coupling the target with opposite spins to the same hnal 
state. The interference between these amplitudes leads to the SSA. Final-state interaction 
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(FSI) from the gluon exchange between the outgoing quark and the target spectator system 
has been used to explain the observed results for the azimuthal SSAs of the proton. In 
particular, we have studied the dependence of transverse azimuthal spin asymmetry {Vy) 
as well as the longitudinal asymmetry (KVy) on light cone momentum fraction A and 
momentum carried by the outgoing quark r^. The implications of such a model have been 
investigated in detail by considering hxed and running coupling constants. It is found that 
the transverse as well as the longitudinal asymmetries change sign at A ~ 0.3 as the value 
of the A increases which may be possibly due to the qq pairs being readily produced at 
this momentum fraction which give rise to sea quarks affecting the sign of asymmetry at 
this point. Several experiments are contemplating the possibility of performing the high 
precision measurements over a wide x region in the near future. Further, Vy and KVy do 
not change sign when we £x A and vary the quark transverse momentum of the outgoing 
quark relative to photon directionr_i_ and decreases as ai A r^. The approach presented 
here can also be applicable to other hadronic inclusive reactions. 

The FSIs also produce a complex phase which can be included in the LFWFs to calculate 
the Sivers and Boer-Mulders distribution functions of the nucleon. We have also computed 
the Sivers and Boer-Mulders distribution function by inducing the spin-dependent complex 
phases to the LFWFs. Sivers function is given by the overlap of the wave functions having 
opposite proton spin states but the same quark spin state whereas the Boer-Mulders function 
has the same proton spin state but opposite quark spin state. We find an overall negative 
sign for both Sivers and Boer-Mulders distribution functions when x ^ 0.3. However, the 
signs reverse for x ^ 0.3. The possibility of a node is also in agreement with a few other 
model calculations and can be related to Qiu-Sterman function. 
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